Persistence in Cluster— Cluster Aggregation 
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Persistence is considered in diffusion-limited cluster-cluster aggregation, in one dimension and 
when the diffusion coefficient of a cluster depends on its size s as D(s) ~ s 7 . The empty and filled 
site persistences are defined as the probabilities, that a site has been either empty or covered by 
a cluster all the time whereas the cluster persistence gives the probability of a cluster to remain 
intact. The filled site one is nonuniversal. The empty site and cluster persistences are found to 
be universal, as supported by analytical arguments and simulations. The empty site case decays 
algebraically with the exponent 9e = 2/(2 — 7). The cluster persistence is related to the small s 
behavior of the cluster size distribution and behaves also algebraically for < 7 < 2 while for 7 < 
the behavior is stretched exponential. In the scaling limit t — > 00 and K(t) — > 00 with t/K(t) fixed 
the distribution of intervals of size k between persistent regions scales as n(k;t) = K~ 2 f(k/K), 
where K(t) ~ t is the average interval size and f(y) = e~ y . For finite t the scaling is poor for 
< t z , due to the insufficient separation of the two length scales: the distances between clusters, t z , 
and that between persistent regions, t e . For the size distribution of persistent regions the time and 
size dependences separate, the latter being independent of the diffusion exponent 7 but depending 
on the initial cluster size distribution. 

PACS numbers: 05.40.-a, 82.20.Mj, 05.50.+q, 05.70.Ln, 02.50.Ey 



I. INTRODUCTION 

Persistence in dynamical systems has attained con- 
siderable interest recently jy. It was originally intro- 
duced for a simple diffusion process ]^|| but since then 
it has been studied, for example, in spin systems [i^l2[, 
reaction-diffusion systems Jl3| -[L9[ , voter model [20pl| 
and for interfaces |22]-|25| ] . It has also been measured 
experimentally pfj|-[30|] , and a few exact results exist 

The persistence is usually defined as the probability 
P(t) that at a fixed point in space a fluctuating nonequi- 
librium field 4>{x; t) does not change sign upto time t, 
i.e., the probability that sgn[<fi(x;t) — (4>(x;t))] remains 
unchanged. In many systems the probability decays alge- 
braically, P(t) ~ t~ e , with a nontrivial persistence expo- 
nent 9. The significance of the phenomenon stems from 
the fact that the exponent is not, in general, related to 
the usual static or dynamic exponents. This this in turn 
implies that not necessarily all of the properties of the 
system are characterized by a single length scale. 

In this article we study various definitions of persis- 
tence in an aggregation process, in the disguise of the 
one-dimensional diffusion-limited cluster-cluster aggre- 
gation model (DLCA) j35|, with each cluster diffusing 
with a size dependent diffusion coefficient, D(s) ~ s 7 
with 7 < 2. This model is in many respects suitable 
for studying persistence as it possesses a rich dynamical 
behavior but at the same time it is simple enough to of- 
ten allow for an analytical treatment [^6|. Moreover in 
DLCA one meets immediately the possibility of defining 
several persistent quantities, each of which describes a 
different aspect of the aggregation process. 

The persistence probabilities considered in this work 



are i) the probability of a cluster to remain unaggregated 
(cluster persistence), ii) the probability that a site has 
been empty (empty site persistence) and iii) filled (filled 
site persistence) upto time t. Notice that all these are 
local quantities. It is possible, in analogy to e.g. the 
contact process j37j, to work starting from definitions 
that involve a global quantity like the average cluster 
size. Thus one could compare the ensemble average of 
the cluster size to that in any particular sample and in- 
vestigate the sign changes of their difference. Likewise 
one could attempt the same, a comparison to the average 
cluster size, for all individual clusters p8|. These would 
be, in practice, much harder to study numerically than 
i) to iii) above which also have the pleasant aspect of be- 
ing, possibly, experimentally relevant. The comparisons 
between the three quantities points out the underlying 
mechanisms, which make some persistent quantities to 
be universal while the others are not. 

Recently Manoj and Ray claimed that the dynamic 
exponent associated with the growth of the intervals 
between two consecutive persistent sites in a diffusion- 
annihilation model A + A — > shows nonuniversal behav- 
ior, i.e., the exponent depends on the initial concentra- 
tion of particles This started a discussion |]l5| , |l6| , 
concerning for example whether it matters if only the 
empty sites are taken to be persistent at t — instead of 
the entire lattice. The problem of universality remained 
unsolved. Later on, the same authors argued that the 
reason might be the competition between the two length 
scales, the diffusive scale Co and the persistence length 
scale C p |jl7|,[l8). Very recently, Bray and O'Donoghue 
studied a similar problem in the g-state Potts model [[| . 
They found no evidence on the concentration depen- 
dence and showed, that the characteristic length scale 
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C c controlling the distribution of the intervals between 
persistent sites, is given by the maximum of the two: 
C c = max{£fl,£j,}. Furthermore they argued that the 
poor scaling for small values of q is due to the insuffi- 
cient separation of these two length scales. These issues 
are easier to resolve in the DLCA where the length scale 
separation is more pronounced. 

First we discuss the three persistence definitions with 
size independent diffusion coefficients (7 = 0) in order to 
clarify the universality of these quantities. The filled site 
persistence turns out to be nonuniversal whereas the two 
others decay similarly independent of the concentration 
and initial conditions. This difference between empty 
and filled site persistence is in clear contrast to all the 
previous persistence studies of, for example, coarsening 
in the Ising model 0, where the symmetry between up 
and down spins is not broken. 

Thereafter we concentrate on the two universal per- 
sistent quantities and consider the influence of size de- 
pendent diffusion (7 7^ 0). We discover that it is pos- 
sible to have situations, in which the other persistence 
decays faster than any power-law in time whereas the 
other shows algebraic behavior. At the end we study the 
distributions of the persistent regions and the intervals 
between these. The region size distributions depend on 
the initial cluster size distribution. The interval size dis- 
tributions are universal, simple exponentials, and their 
scaling nicely demonstrates the effect of the presence of 
the two length scales, C p and Cjj. 

This paper starts by introducing the model and de- 
scribing the quantities of interest in section ||. The 
known results for the persistence probabilities are pre- 
sented in the beginning of section [II. In the rest of 
that section each persistence probability is analyzed sep- 
arately. Section IV discusses the scaling of the region and 



interval size distributions. The dependence of the persis- 
tent quantities on concentration and initial conditions are 
studied by simulations in the beginning of section 0. The 
end of that section shows the numerical results for the re- 
gion and interval size distributions. Section |Vl] concludes 
the paper. 
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connected via nearest neighbor occupancy belong to the 
same cluster. Each cluster performs a random walk such 
that the diffusion coefficient of a cluster of size s is given 
by D(s) = Dis 1 . The positive constant D\ sets the time 
scale and it is irrelevant regarding the dynamic scaling 
properties of the system which are controlled by the 
diffusion exponent 7 ]39|-|42]|. Here we concentrate on 
7 < 2 for which the growth of the average cluster size is 
algebraic. For 7 > 2 the system will gel, i.e., there will 
be an infinite cluster at a finite time. When two clusters 
collide, they aggregate irreversibly together. If the sizes 
of the colliding clusters were s and s', the aggregate is of 
size s + s' and its diffusion constant D(s + s'). Figure |l| 
visualizes the model. 



The one-dimensional DLCA is closely related to other 
models when the diffusion coefficient of a cluster is in- 
dependent of its size, i.e., 7 = 0. For example, con- 
sider the zero temperature one-dimensional g-state Potts 
model with Glauber dynamics. The different spin species 
are separated by domain walls, each of which performs a 
simple random walk. By representing each such a walker 
by a particle A, the two particles colliding react according 
to 



A + A 



A with probability 
with probability 



(g-2)/(«-i) 

1/(9-1)- 



(1) 



FIG. 1. Illustration of the one-dimensional DLCA model. 



In the limit q — ► 00 this reduces to the reaction-diffusion 
model A + A ^ A, which can be solved exactly pq |. The 
same is true for particles with a finite mass Ai + Aj — > 
Ai+j @|. The only difference between this and the 
DLCA model is the finite extent of clusters in the DLCA. 
It has no relevance on the cluster size distribution but 
the difference may be essential for persistence. This is 
to be expected since i) the g-state Potts model can be 
mapped to the Ising model (q = 2) with an initial frac- 
tion 1/q of up spins || and ii) the persistence exponent 
of the Ising model is known to depend on the concentra- 
tion [HQ]. Note, that in the Potts model the usual defi- 
nition for persistence characterizes the fraction of spins, 
which have remained in their original state. 



II. MODEL AND QUANTITIES OF INTEREST 

The diffusion-limited cluster-cluster aggregation 
model is here considered on a one-dimensional lattice 
of L sites with periodic boundary conditions. Initially 
the lattice is filled upto a concentration 4> such that 
occupied lattice sites correspond to particles and sites 



There are many ways to define persistence in cluster- 
cluster aggregation. Perhaps the most closest to the def- 
inition used in the context of Ising systems would be the 
one which relates to the question of site persistence: what 
is the probability, that a site has been either covered or 
uncovered by a cluster for all times [0,i]? However, in 
DLCA there are also other natural candidates. 
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FIG. 2. A space-time plot of the dynamics in DLCA for 
7 = 0. The clusters are denoted by light gray, the persistent 
clusters by black and the persistent empty regions by dark 
grey. 

The persistence probabilities studied in this work are 
the following: 



• Cluster persistence: 
has not aggregated, 



the probability that a cluster 

P C (t)~t- 6 c. 



• Filled site persistence: the probability that a site 
originally covered by a cluster has been covered by 
it all the time, P F {t) ~ t~ 9F . 

• Empty site persistence: the probability that an 
originally empty site has never been occupied by 
a cluster, P E (t) ~ t~ 0E . 

When the persistence probabilities decay algebraically, 
the associated decay exponents 6c, Of, and E become 
persistence exponents. Figure ^ shows an example of the 
dynamics together with the persistent clusters and the 
persistent empty sites for a small system. 

In the following we will concentrate only on the three 
definitions given above. However, some other persis- 
tences can be directly obtained from the results given. 
For example, the aforementioned site persistence prob- 
ability is defined as the probability that a site has re- 
mained in the same state, i.e. either filled or empty, all 
the time and is given by P s (t) = Pp(t) + P E (t). The 
survival of domains follows from the studies of the per- 
sistent region size distribution (defined below) . This is a 



special case of the so called adaptive persistence 44 



Furthermore one could consider the particle persis- 
tence, which is the probability that a site has remained 
either empty or covered by the same particle all the time. 
Recall that a cluster of size s consist of s particles. The 
probability that a cluster has not moved before time t 
decays exponentially making particle persistence a rather 
trivial quantity. 



The cluster persistence differs drastically from the 
other persistence definitions. It is not a quantity de- 
fined per a fixed site on the lattice as the others are 
but is a property associated with each cluster. Due to 
that also its behavior is rather different than that of the 
other persistences considered. However, it is a rather 
similar, though not exactly the same, quantity as the 
domain persistence in the Ising model H0. For point- 
like clusters it leads to a study of survival of coalesc- 
ing random walks. A similar study has very recently 
been carried out for non-interacting walks and for the 
reaction-diffusion system A + B — > in jl9| and for the 
annihilation-coalescence reaction-diffusion system corre- 
sponding to the Potts model [Eq.(j|)] in 

In addition to the persistence probabilities we will also 
study distributions associated with the persistent sites 
and the distances between them. The terminology and 
notation used are as follows. The word region is reserved 
for a bunch of consecutive persistent sites. The distances 
between regions, i.e., between two consecutive persistent 
sites, is called an interval. The word cluster has the obvi- 
ous meaning. The number of clusters of size s per lattice 
site at time t is n s (t) with the normalization ^ s sn s = 1. 
Region size is denoted by I and the number of regions of 
size I (per site) by px(l;t). The subscript is like for the 
persistence probabilities and it refers to the persistence 
definition used: X 6 {C,F,E}. When using the contin- 
uum description we use the symbol r instead of /. The 
letter k labels the interval sizes and the corresponding 
distribution function is nx(k;t). The number densities 
of persistent sites (the persistence probability) and non- 
persistent ones are denoted by capital letters Px(t) and 
Nx(t), respectively. 

As an example, consider intervals between persis- 
tent empty sites and their distribution function n E (k;t). 
The corresponding distribution of persistent regions is 
p E (l;t). The number densities are obtained by summing 



N E (t)=J2 kn E (k;t) 
fc=i 

oo 

P E (t) = J2 lp E (l;t). 



(2) 
(3) 



Obviously these two are related by the equation P E (t) = 
1 — N E (t). Similar formulas apply to other persistence 
definitions, too, except in the cluster persistence case 
equation (||) is replaced by 



Pc(t) = J2p c (l;t) 



(4) 



and naturally P c {t) ^ 1 - N c (t). 



III. PERSISTENCE PROBABILITIES 

In this section we first list the known exact results for 
7 = 0. To discuss the cluster persistence we use the 
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known results of the dynamic scaling of the cluster size 
distribution. We also briefly comment on how the cluster 
persistence exponent is related to the polydispersity ex- 
ponent, which characterizes the behavior of the scaling 
function of the cluster size distribution. Then we ap- 
ply mean-field type continuum methods to analyze the 
empty site persistences for 7^0 where no exact results 
are available. The filled site persistence is analyzed in a 
similar mean- field fashion in the limits and </> — > 1. 

This section ends to a discussion on the concentration de- 
pendence of the filled site persistence. 



A. Known Exact Results 

For the cluster persistence the finite extent of the clus- 
ters does not matter and we may replace them by point- 
like random walkers. Then we can utilize the results ob- 
tained for the the reaction-diffusion model Ai + Aj — > 
Ai+j. For spatially uncorrelated initial mass distribu- 
tion n s (0) = S s s the cluster size distribution is given by 
n s (t) = 2se - *ia (£)/£) where f = 4Di* and I a {£) is the 
modified Bessel function J3f|. As ri\(t) ~ i -3 / 2 the clus- 
ter persistence exponent #c(7 = 0) = 3/2 in this case. 
This result is also valid independent of the initial condi- 
tion as long as the cluster size distribution decays rapidly 
for large sizes. 

The result 9 c (0) = 3/2 can also be obtained without 
solving the cluster size distribution for general initial size 
distribution of clusters. Consider a cluster initially at the 
origin. We want to find the probability that it has not 
collided with neither of its neighboring clusters before 
time t. As the diffusion coefficient of a cluster does not 
depend on its size, we do not have to care if the neigh- 
boring clusters have collided with other clusters or not. 
Therefore the problem reduces to the survival probabil- 
ity of three annihilating random walkers. This can solved 
in several ways (45||| with the result 9c{l = 0) = 3/2. 
Unfortunately, these methods rely on the indistinguisha- 
bility of clusters and can not simply be generalized to 
the general case of size dependent diffusion coefficients. 
In fact, even the survival of three annihilating random 
walkers with constant but non-equal diffusion coefficients 
is highly nontrivial [^6|. In DLCA it is the specific in- 
teraction between clusters, the aggregation at contact, 
that makes the problem simple. The case where a single 
cluster is hunted by independent clusters is in fact more 
complicated than our situation with aggregating clusters 
(see e.g. @). 

Let us now turn to site persistence when 7 = 0. For 
the empty site persistence the finite extent of the clusters 
does not matter since a persistent interval is deleted by 
the edges of neighboring clusters. For finite size clusters 
the empty intervals are shorter but the length distribu- 
tion remains the same. The movement of the edges is 
uncorrelated and unaffected by the collisions with other 
clusters. These facts make the problem exactly solvable 



with the result 9 E (j = 0) = 1 f§. Note, that although 
the movement of edges is uncorrelated for 7 7^ 0, too, 
one can no more ignore the collisions of the neighboring 
clusters with other clusters. This makes the exact calcu- 
lation complicated if not intractable. The same applies 
to the filled site persistence even for 7 = 0. 



B. Cluster persistence 

Before considering cluster persistence we briefly reca- 
pitulate the main results concerning the dynamic scaling 
in the DLCA. Let the number of clusters of size s per lat- 
tice site at time t be n s (t). The cluster size distribution 
scales as 



n s (t) = S(t)- 2 f(s/S(t)), 



(5) 



where S(t) is the average cluster size and f{x) the scal- 
ing function pEfl . The diffusive length scale, Cd, is 
proportional to the average cluster size which grows as 
S(t) ~ t z , where z is the dynamic exponent. For a given 
7 both z and the scaling function are universal. In one di- 
mension the dynamic exponent is given by z = 1/(2 — 7) 
for 7 < 2 and the system gels in a finite time for 
7 > 2 For small argument values [x — > + ) the 

scaling function f(x) decays as x~ T for < 7 < 2 and as 
exp(— aH^I) for 7 < Though universal, both 

the exponents r and fi depend nontrivially on 7. From 
the small x behavior of the scaling function it follows that 
n s (t) ~ s- T t- w forl«s« S{t), where w = (2 - t)z. 

The scaling theory further states that for the DLCA all 
the cluster densities decay in a similar manner at large 
times: n s (t)/ni(t) — > b s as t — > 00 [Esl. Here b s is a 
time independent constant. This result together with 
n s (t) r~j f w can be directly applied to cluster persistence 
with a monodisperse initial condition n s (0) = <5 s . Sq . For 
simplicity, let all the clusters be monomers, i.e., of unit 
size, at time t — 0. Then the cluster persistence proba- 
bility is equal to n\{t) ~ t~ w giving 



9 C - (2 - t)z. 



(6) 



The exact exponents r = — 1, z = 1/2 j36|, and w — 
3/2 |4f],|^] for 7 = satisfy this relation. Simulations 
show that similarly than z and r also the persistence ex- 
ponent is indepe nden t of initial conditions and concen- 
tration (see Sec. VB ). 

The result 6c = (2 — t)z deserves a few comments. 
First, since the cluster persistence is a cluster property, 
it is natural that the persistence exponent is related to 
the cluster size distribution. Since a persistent cluster has 
not collided with others, it will definitely be smaller than 
the average cluster size at late times. Thus the persis- 
tent clusters are the ones represented by the small s-tail 
of the cluster size distribution, which is described by the 
polydispersity exponent r. On the other hand, if there is 
any scaling relation between the cluster persistence expo- 
nent and the traditional exponents, this should include 
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the dynamic exponent. On the basis of this reasoning, 
the above result is reasonable. However, most of the 
scaling results quoted above are mean-field of nature and 
therefore Eq. (g) has to be validated by simulations. The 
same scaling relation has been found also in the diffusion- 
annihilation model discussed in the Introduction [ff8f . 

It follows from the scaling relation that the cluster per- 
sistence exponent is as universal as the exponents z, w, 
and r. However, it does not result in a value for the 
persistence exponent since r is a priori unknown. Calcu- 
lating it for DLCA and other similar processes has turned 
out to be a challenging problem since it is given in terms 
of integrals over the scaling function f(x) itself p9| , p9[ . 
Also determining the value of r reliably in simulations is 
hard. On the other hand, if one could calculate the per- 
sistence exponent and validate the scaling relation one 
would also obtain r. 

In fact, the persistence exponent can be calculated in 
a mean-field fashion by replacing the clusters surround- 
ing a persistent cluster by random walkers which diffuse 
as an average cluster would diffuse. A similar approach 
works ex treme ly well also for the empty site persistence 
(see Sec. IIIC). Here it results in studying the survival 
of three random walkers, whose diffusion coefficients are 
time-dependent. The meaning of a time-dependent dif- 
fusion coefficient in this context is that a random walker 
will diffuse with a constant diffusion coefficient in another 
time scale. The detailed analysis and the numerical ver- 
ification of Eq. (H) will be presented elsewhere j5Tj] and 
we just give the result here. 

It turns out that the cluster persistence probability 



Pc(t) 




(7) 



with C > a constant. j3s fits well to the expression 
f3s = 2/3(1 — 2z). Here the completely different behavior 
of the cluster persistence for 7 < and 7 > is similar to 
that of the cluster size distribution which changes char- 
acter at 7 = as noted in the beginning of this Section. 
The persistence exponent is discontinuous as 7 — * + . 
Quite unexpectedly, since 9 c (0 < 7 < |) < #c(0), it 
is easier for a cluster to remain persistent if the other 
clusters will diffuse faster! This somewhat paradoxical 
result can be understood on the basis that for a random 
walker, which is hunted by other random walkers, the op- 
timal strategy to survive is to remain stationary (see |5lj] 
for a similar situation in the case of noninteracting walk- 
ers hunting a single random walker) . In the limit 7 — * 00 
the stretching exponent [3s — ► 2/3 is also "discontinu- 
ous" , since the persistence for a particle bounded by two 
immobile neighbors decays simply exponentially. 

The result (0) together with Eq. (||) implies that also 
the polydispersity exponent is discontinuous at 7 = 
and, more interestingly, that 7- (7) = for < 7 < 2. Al- 
though this kind of a "superuniversality" is remarkable 
in its own right, it is more important, that by studying 



the cluster persistence one learns something about the 
polydispersity exponent, i.e., the scaling function. This 
connection may exist in other aggregation models offer- 
ing a new way to approach the problem of computing the 
small size tail of the scaling function through a perhaps 
simpler quantity, the cluster persistence. 



C. Empty site persistence 



Next we give a heuristic argument for obtaining the 
empty site persistence exponent for arbitrary 7. Since 
the clusters on both sides of a persistent empty region 
are independent, we are led to consider the maximum 
excursion of a single random walker. The only complica- 
tion comes from the fact that for 7 7^ the collisions of 
a cluster with other clusters will change its diffusivity. 

Instead of the real process we consider each cluster 
to have a time- dependent diffusion coefficient D(t) — 
D t^ z ~ S(t)f with z = 1/(2-7). Then the proba- 
bility, P(x;t), to find a cluster at position x at time t 
obeys a diffusion equation 



d t P(x;t) = D t Jz d*P(x;t). 



A time transformation 



T (t) = J^—t^ z+1 
jz + 1 



(8) 



(9) 



reduces equation (||) to an ordinary diffusion equation 
with the diffusion constant D = 1. The persistence of 
the empty space between particles diffusing with a con- 
stant D has recently been considered in ||] and we just 
quote the main results here. 

In the long time limit the probability density, 
PE(r;t\s), that a persistent empty region (originally of 
size s) has size r at time T is given by 



p E (r;T\s) 



s — r 



(10) 



and the probability that a cluster survives upto time t is 



P E {t\s) = 



2ttT ttDq t 2z 



(11) 



For a general initial length distribution of clusters n s (0) 
the result will remain the same except s 2 will be replaced 
by the average over the initial length distribution (s 2 ). 

These considerations show that the persistence expo- 
nent 9e{i) = 2z = 2/(2 — 7). This agrees with the exact 
result 9e(0) = 1. Furthermore the result is indepen- 
dent of the initial spatial distribution or concentration. 
All these implications based on the simp l ified process are 
confirmed by simulations (see Sees. V A- VC). 
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FIG. 3. A space-time plot for concentration (j> = 0.8 and 
7 = 0. Cluster are denoted by grey and the persistent empty 
sites by black. 



D. Filled site persistence 

The behavior of the filled site persistence is different 
at low and high concentrations. At low concentration 
a cluster contains typically only one persistent region 
which is usually destroyed before the cluster aggregates 
with a neighbor. At high concentrations a large cluster 
has many persistent regions of various sizes inside it due 
to the aggregation of clusters (Fig. ||). We now analyze 
these two limits separately. 

In the low concentration limit we consider the persis- 
tence problem in continuum, for simplicity. This should 
be valid for clusters of initial length Iq 3> 1. For 
the time required for a cluster to move its own length is 
much smaller than the time required to overcome the dis- 
tance between clusters. Therefore one could assume that 
at low concentrations collisions between clusters do not 



matter. We will see in section V A that this is true only 
upto some crossover-time. This time obviously diverges 
in the limit <f> — > 0. However, to obtain more insight to 
the empty site persistence we will first ignore the colli- 
sions in the following analysis. 



Xmin 


Xmax R(t) 


^min 


^max 



















nonpersistent sites 

FIG. 4. The length of the nonpersistent filled part is equal 
to the span R(t) = x m in(t) + £ max (f) of the middle point of 
a cluster. The length of the persistent region is so — R(t). 
The cluster position at t = is denoted by a thick rectangle, 
whose middle position is marked by a dotted line. The dashed 
lines show the maximum excursions of the cluster. 



When clusters do not collide, the persistent sites under 
different clusters are destroyed by independent processes. 
A single diffusing cluster destroys persistent area at both 
of its edges. As figure || shows, the length of the nonper- 
sistent filled part inside an unaggregated cluster is equal 
to the span of the middle point. The span, R(t), of a ran- 
dom walk is defined as R(t) = x m - ln (t) + x max (i), where 
Xmin(t) and x max (t) are the maximum displacements of 
the walk in the negative and positive directions at time 
i, respectively. 

The probability distribution for the span of an unlim- 
ited random walk is given by |53] 



w(R;t) = 



£(-iy+Vexp(- 



j 2 R 2 

2Dt 



(12) 



In our case the maximum span is limited by the size of 
the cluster, s, and the probability distribution p F (r; t\s) 
that an interval of length r of a cluster initially of size s 
has survived upto time t is 



p F (r;t\s) 



V2wDt 



£(-iy+y ex p(- 



j 2 (s - rf 



2Dt 



(13) 



for r < s and zero otherwise. 

We are interested in the asymptotic long time behav- 
ior and rather than working directly with the formula in 
Eq. jl3|) it is useful to transform it to a more tractable 
form as follows. Writing the sum in Eq. (13) in the form 
Ej*l-oo and applying the Poisson sum formula |53"f 



m— — oo 



E g*= E 

leads to 

p F (r;t\s) 



8Dt 



(s — r) 3 
x exp 



E 



da; G{x)e 



TT 2 Dt(2m+l) 2 



(14) 



(s - r) 2 
m=0 L v ; 

ir 2 Dt(2m+ l) 2 

2(s-r) 2 



(15) 



In Eq. ( |l4| ) G(x) is an analytic function having the same 
values at all integers as Gj. 

From Eq. (|l5|) it is easy to show that the probability 
of finding persistent sites inside the cluster decays expo- 
nentially at large times: 



p F (t\s) 



dr p F (r;t\s) 



E 

m=0 

x exp 
8Dt 



[s 2 +ir 2 Dt(2m+ 1)] 



s 2 TT 2 (2m + l) 2 

Tr 2 Dt(2m+ l) 2 
2? 
n 2 Dt 



■ exp 



2s 2 



(16) 
(17) 



G 



The persistence probability is obtained by integration 

p F (t) = / ds PF (t\s)n s (0), (18) 



so that, in general, the decay will depend on n s (0). 

For n s (0) — S(s — Sq) it is given by an exponential ex- 
pression [Eq. (|l7|)l but the change in the decay can be 
shown using, for example, the initial distribution 



2 ( s 
n s (Q) = exp 2 

7TS \ TTSq 



(19) 



whose mean and standard deviation are (s) = sq and 



(T s = yj (s 2 ) — (s) 2 — soy^jr — 2)/2, respectively. Using 
the approximation (|17|) we get a stretched exponential 



PF(t) 



7T 3 S 2 



For the Poisson initial distribution, n s (0) = s e 
we can apply the saddle point method to give 

p F (t) : 



'128Dt / 3 /n 2 Dt\ 

2~ ex P 



Sirs 



(i 



1/3N 



2 v 4 J 



(20) 



s/s 



(21) 



Both of these examples show a stretched exponential 
decay for the persistence probability, which we believe 
to be true quite generally. The stretching exponent de- 
pends on the initial condition and is therefore nonuni- 
versal. The long time behavior, in the approximation 
ignoring collisions, is governed by the ratio F)tj s\ and 
we also immediately see that the decay depends on con- 
centration due to the factor sq in the exponential. 

The approximative span distribution can be also used 
to calculate the mean size of persistent regions. First 
write 



(r(t\s)) = (s-R(t\ S )) = s-(R(t\s)) 



(22) 



and use the normalized approximate form for the span 
distribution 



'S(R-s) 1- (l + £s- 2 )e-«/ s 
w N (R;t\s) = { + 2 £ 2 i?- 5 e-«/ fl2 
.0 



;R< s 

;R>s, 
(23) 



where £ = ir 2 Dt/2, to calculate (R(t\s)) from which 



(r(t\s)} « s exp < — 



TT 2 Dt 

2s 2 



(24) 



For the simplest case of a fixed initial size so [n s (0) 
5{s — sq)] the mean length of surviving regions is 



PF(t\s ) J 



dr rp F (r;t\s ) 



2*g 
ir 2 Dt 



(25) 



This is just an example and naturally the decay depends 
on the initial distribution n s (0). 

In the high concentration limit, <f> — * 1, we adopt an- 
other mean-field type approach: we consider a deter- 
ministic model combined with scaling arguments. Let 
initially the average cluster and empty interval size be 
So and size do, respectively. The concentration (f> — 
so/(«o + do) ~ 1 f° r «0 > do- 

Now consider the doubling times tx, t%, . . . , ti, . . . at 
which the average cluster and interval sizes are Si = 2 1 sq 
and di — 2 l do, respectively. At each step n the dou- 
bled cluster is constructed as follows. First, d„_i/2 sites 
(these does not have to persistent but they may be) from 
both ends of the cluster are made nonpersistent. Sec- 
ond, the resulting cluster is duplicated. This procedure 
is repeated. It is quite easy to convince oneself that the 
probability of finding persistent sites, p(n), at step n de- 
cays like p(n) - e - Q W™ for high enough values of 6. 
Since s„ = 2"s ~ t z , it follows that p(t) - ^ 

According to this simple and crude argument the filled 
site persistence probability decays algebraically for large 
enough concentrations. The persistent sites are sweeped 
by the domain walls (clusters' edges), which get destroyed 
at aggregation. Since the probability to be touched by 
a domain wall depends naturally on the density of the 
walls in the system, the persistence exponent depends 
on concentration implying nonuniversality. Simulations 
qualitatively agree with this behavior. We have not made 
any quantitative comparison since in this paper we con- 
centrate more on the universal properties of the cluster- 
cluster aggregation and since being non-universal the de- 
cay of site persistence is harder to compute analytically. 
The simulations furthermore show that the persistence 
probability decays algebraically for low concentrations, 
too. The reason for the deviation from the span argu- 
ment lies in the approximation, which neglects the colli- 
sions between clusters. We return to this issue in a more 



detail in the section V A 



IV. 



PERSISTENT INTERVAL AND REGION 
DISTRIBUTIONS 



Having analyzed the various persistence quantities we 
now turn to the distributions of persistent regions and 
the intervals between them. First we consider the in- 
terval size distributions between consecutive persistent 
sites and perform a scaling analysis. A natural scaling 
assumption for these distributions is 



n x (k;t) = K x (t)~ a f x (k/K x (t)), 



(26) 



where Kx(t) ~ t 8x is the average interval size and X 
denotes any of the choices C, E, or F. 
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Inserting the scaling form (|26|) to the equation (|2j) and 
replacing the sum by an integral (valid in the long time 
limit) results in 



V. SIMULATIONS 



N x {t)=K x (t) 



2-a 



yfx(v)- 



(27) 



As Px{t) — * for t — > oo it follows from the relation 
P x {t) = 1 - N x {t) that iVje(i) — ► 1 as t — ► oo. The only 
way to keep the right hand side of Eq. ( p?j) constant is 
to require a = 2, as a direct consequence of mass conser- 
vation. Note, that the above argument does not require 
the persistence probability Px(t) to decay algebraically. 

Since the persistence exponents are larger than the dy- 
namic exponent, 6x > z, the persistence length scale will 
be much larger than the diffusive one at large times. This 
means that the persistent regions are well separated and 
that they are destroyed by uncorrelated processes since 
the correlations grow only as t z . Each persistent region 
survives independently of the others and in the scaling 
limit the scaling functions are therefore simple exponen- 
tials 



n x (k;t) = K x (t) 



-2 e -k/K x {t)_ 



(28) 



Next let us turn to the size distribution of persistent 
regions. As the filled site persistence shows nonuniversal 
behavior, we will here concentrate only on the distribu- 
tions concerning the empty site and cluster persistence. 
The former has already been analyzed in section III . For 



example, equation ( |10|) gives the persistent empty region 
distribution at large times for a monodisperse initial dis- 
tribution and for other initial distributions n s (0) it is 
obtained as 



Pe(t;T) 



ds p E (r; T|s)ra a (0) = -^e 



-r/so 



(29) 



where T is given by equation (|9J) and the last form cor- 
responds to the initial distribution ri s (0) = s^e^ 8 ' 80 . 
The dependence on the diffusion exponent enters only 
through the time scale T ~ t 2 /^ 2 ~ 1 \ 

More important is that the spatial and time depen- 
dence in p E (r;t) are decoupled. From this it follows 
that the average size of the persistent empty regions, 
L E (t), is a constant at large times. For the monodis- 
perse and Poisson initial conditions we get L™ ono (i) = 
Pe^so)- 1 f°° dr r PE (r;T\so) = so/3 and L r E and (i) = 
/ °°dr rp E (r;T)/ J °° dr p E (r\T) = s , respectively. 
These b oth a re independent of time. The simulations 
(see Sec. VE ) confirm this. 

Unfortunately we have not been able to similarly cal- 
culate the size distribution of persistent clusters, pc{l\ t). 
However, for initial cluster size distribution n s (0) = <5 So ,s 
this distribution remains unaltered. The same is true also 
for size independent diffusion coefficients (7 = 0) no mat- 
ter what the initial distribution is. This does not have 
to be true when 7 7^ 0. For instance, there may be non- 
trivial correlations between a persistent cluster and the 
sizes of clusters surroundin g it. The scaling of pc(l',t) is 
briefly discussed in section | 



The simulations are done on a one-dimensional lattice 
with periodic boundary conditions. In all the simula- 
tions D\ = 1 fixing the time scale. The system sizes 
range from 5 ■ I0 4 to 1.5 • 10 6 , the data are averaged 
over 1000 — 50000 realizations. The algorithm used is as 
follows. A cluster is selected randomly and the time is 
incremented by N^)^ 1 D^^, where N(t) is the number 
of clusters at time t and -D max is the maximum diffu- 
sivity of any of the clusters in the system at that time. 
If x < D(s)/ Z? ma x, where x is an uniformly distributed 
random number in the interval [0, 1], the cluster is moved 
with equal probability either to the left or to the right. 
Otherwise a new cluster is selected and the process is re- 
peated. If a moved cluster collides with another one, the 
two clusters are aggregated together. 



We test the sensitivity of the per sisten ce probabilities 
against concentration chan ges ( Sec. V A) and two differ- 
ent initial conditions (Sec. |VB| ). The first initial condi- 
tion used is random, which means that each lattice site is 
filled with probability p. The other one is deterministic 
and monodisperse such that clusters of a given size Sq 
are separated form neighboring clusters by a distance Iq 
resulting in the concentration <j> = so/(sq + Iq). 



After discussing the universality or nonuniversality of 
the persistences we concentrate on the dependence of the 
persistence exponents on the diffusion exponent 7 in sec- 
tion |V C| . Ther eafte r we consider the size interval dis- 
tributions (Sec. |VD| ). We demonstrate how the poor 
dynamic scaling of the interval distributions is a conse- 
quence of the two competing length scales C E and C p . 
The scaling function of the region size distributio ns is 
shown to be of the simple exponential form in section V E . 
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FIG. 5. The persistence probabilities a) Pe(£), b) Pp{t), 
and c) Pc{i) for size independent diffusion coefficients (7 = 0) 
and for concentrations (p = 0.20 (O), 0.50 (□), and 0.80 (V). 
Results from simulations using random and monodisperse 
(so = I) initial conditions are denoted by open and filled 
symbols, respectively. Pluses (+) are obtained with monodis- 
perse initial conditions for <j> = 0.20 and so = 10. The insets 
show the running exponents. 



A. Concentration dependence 

We first analyze which persistence probabilities depend 
on concentration. We will present data only for size 
independent diffusion coefficients but we have checked 
that for 7/O the persistence probabilities behave in 
a similar way. Figure |^ shows the cluster, empty site, 
and filled site persistence probabilities for three differ- 
ent concentrations. A change in the concentration af- 
fects only the amplitudes of the cluster and empty site 
persistence distributions. The numerical estimates for 
the concentration independent persistence exponents ob- 
tained from the saturated part of the running exponents 
are C (O) = 1.48 ± 0.03 and 9 E {0) = 1.00 ± 0.02. These 
both are in an excellent agreement with the exact results 

M 



; (0) = § and 9 E (0) = 1. 
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FIG. 6. a) Filled site persistence probability and b) aver- 
age cluster size of those clusters which contain persistent sites 
for deterministic initial conditions with 7 = 0, so = 11, and 
4> — 0.2. The solid line is given by equation ( p^ ) and the 
dashed line shows t 1 ^ 2 behavior. 

For the filled site persistence the effect of a concentra- 
tion change is much more pronounced. The distribution 
shows a transition from an algebraic decay to a relatively 
faster one when decreasing the concentration. This is in 
a qualitative agreement with the analysis of section III D 
but the agreement is not quantitative. In fact, the non- 
algebraic decay seen in simulations is only a cross-over 
behavior and eventually the persistence probability will 
follow a power-law. 

The discrepancy between the analysis of Sec. [II D| and 
simulations can be traced to the assumption, that the 



clusters containing persistent filled sites have not aggre- 
gated. This is true only for times smaller than the av- 
erage collision time between clusters, t co \\ ~ l/(D(f) 2 ), 
which indeed diverges for <j> — ► 0. However, for any 
nonzero concentration there will be a time after which 
the collisions become important. In fact, the large time 
persistence probability of filled sites is dominated by the 
clusters, which have collided with other ones. This is il- 
lustrated in Figure ^, which shows both the persistence 
probability and the average size of those clusters which 
contain persistent sites. There is a clear cros sover from 
the behavior given by the analysis of Sec. HID at t ss 300 
to the one for which the collisions are significant. Af- 
ter this crossover time, the clusters including persistent 
empty sites grow similarly as the other ones. This is il- 
lustrated in figure^ b), which shows that the average size 
of these clusters grows like t 1 ^ 2 at late times. 
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FIG. 8. Comparison of the numerically obtained empty 
site persistence exponent 0e (□) to the mean-field result 
6e = 2z = 2/(2 — 7) (solid line) as a function of the diffusion 
exponent 7. 



C. Dependence on 7 



B. Dependence on initial conditions 

Figure ^| also compares the persistence probabilities for 
the two initial conditions: random and monodisperse de- 
terministic. A change in the initial condition does not 
have a significant effect on any of the persistence prob- 
abilities. It only affects the overall amplitude but leaves 
the qualitative behavior unaltered. 

Simulations can of course not prove the universality of 
empty or cluster persistence but together with the ar- 
guments given in section III they strongly support this. 



However, they indisputably show that the filled site per- 
sistence is nonuniversal. Although non-universal, the 
empty site persistence is closely related to the persistence 
of bubbles in soap froths, which show similar concentra- 
tion dependent behavior [ p7| . It would be interesting to 
study the similarities between these two persistences, but 
we will omit it here. 
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FIG. 7. The cluster and empty site persistence probabil- 
ities Pc{i) (filled symbols) and -Pfi(t) (empty symbols) for 
7 = -1.0 (V),0.0 (□), and 0.75 (Q)- 



Thus far we have considered only systems with 7 = 0. 
All the dynamic scaling properties depend on 7 and the 
same is true for the persistence probabilities as figure |?j 
shows. The empty site persistence probability decays 
algebraically for all values of 7 < 2, above which the 
system gels in a finite time. The analytical prediction, 
0e = 2z, obtained by replacing the discrete changes in 
the sizes of clusters besides the persistent one by a diffu- 
sion process, is compared to simulations in figure ||. The 
agreement is excellent, which basically means, that the 
system self-averages away the discreteness of the collision 
times. It also tells that there are no nontrivial correla- 
tions such that the clusters surrounding persistent empty 
sites would grow differently than the other ones. 

The cluster persistence decays algebraically only for 
7 > and faster than any power of t for 7 < 0. The 
nonalgebraic behavior is not surpris ing o n the basis of 
the discussion presented in section IIIB. Although a 



detailed study of the cluste r pers istence confirming the 
results presented in section IIIB will be presented else- 
where []50f| , we briefly comment on the simulational prob- 
lems one can encounter. 

First, as the persistence exponent is related to the 
polydispersity exponent through 9c — (2 — t)z, we are 
faced with same simulation problems than in the stud- 
ies of r (see e.g. 40,4J|). Also the crossover time to the 
asymptotic behavior is large for 7 positive but close to 
zero. In order to see the real asymptotic behavior there 
must be very large clusters in the system. Consider, for 
example, 7 = 0.2, which still is not very small. The diffu- 
sion coefficient of a persistent cluster is of order D\ and 
that of a typical cluster DiS(t)' y . The ratio of these two 
to be of order 10, say, requires S(t) w IO 1 / 7 w 10 5 , which 
corresponds to a crossover time t CTOSS ~ 10 9 . Even reach- 
ing this in simulations would be hard. This problem can, 
of course, be overcome by using a more sophisticated al- 
gorithm. We refer to |5(| for more details on the cluster 
persistence. 
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D. Intervals between Consecutive Persistent Sites 



In section IV we argued that the size distribution of in- 
tervals between persistent quantities would scale accord- 
ing to nx(k;t) = Kx(t)~ 2 fx{k/ 'Kx(t)) with a simple 
exponential scaling function [Eq. (0)]. When the cor- 
responding persistence probability decays algebraically, 
Kx(t) ~ Axt 6x , with a nonuniversal amplitude Ax{4>), 
this can be presented also as nx(k;t) = t~ 26x fx(k/t Bx ). 

The difference between these two scaling forms is that 
for the latter the scaling functions will not overlap each 
other for different concentrations due to the nonuniversal 
amplitude dependence. Therefore we show for clarity the 
scaling plots using the latter formulation. Furthermore, 
we prefer to show the scaling of the complement of the 
cumulative distribution 




FIG. 10. Scaling plot for the distribution of inter- 
vals between two consecutive persistent empty sites at 
t = 2 6 , . . . , 2 11 for <j> = 0.2 and 7 = 0. 



I x (k;t) = ^2<nx(i;t), 

i>k 



(30) 



which is much smoother due to the summation. It is easy 
to see from the scaling of nx(k; t) and Eq. (E8|) that this 
distribution should scale as Ix(k;t) = t~^gx(k/t 9x ) 

1 - — 1 

with g x (x) = A x 



Figure |^ shows that the scaling works for the empty 
site persistence and that the scaling function is an expo- 
nential one indeed. The plots for the other persistence 
definitions are similar and are not shown. The scaling 
function is universal and all the curves in figure ^ would 
overlap each other if one would plot KE(t)I(k;t) as a 
function of k/KE(t). Note, that the diffusion exponent 
7 has no influence on the scaling function. 
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FIG. 9. Scaling plot for the complement of the cumulative 
distribution of intervals between two consecutive persistent 
empty sites at t = 2 9 ,... , 2 14 for concentrations and diffu- 
sion exponents shown in the figure. The value used for the 
persistence exponent Be = 2/(2 — 7). 



Although the summation in Eq. (^0|) smoothens the 
data it at the same time loses information about the 
small k/t® B behavior. This is illustrated in figure [l(^ 
where no summation is done. For k <C t 0E the scaling 
does not work. The reason is the following. The scal- 
ing should work in the limit k — > 00 and t — > 00 with 
y = k/t 0E fixed. Especially, the condition k 3> t z should 
be satisfied for the two length scales of the problem to 
be well separated. Define now a time-dependent y c (t) so 
that the scaling works for y > y c (t). This quantity gets 
smaller at the same rate as the curves in figure [Io| shift to- 
wards zero. Our estimate from the numerical data gives 
y c (t) ~ t"°- 50±0 03 , which is consistent with y c (t) ~ t~ z . 
Thus the poor scaling in figure [LCj for k < t 0B is just a 
manifestation of the finite time behavior with two com- 
peting length scales, Cd and C p . This effect vanishes in 
the scaling limit. A similar, although not as clear, viola- 
tion of scaling induced by the same reason is seen in the 
study of persistence in the g-state Potts model j| . 



11 




FIG. 11. The size distributions of persistent empty regions 
for 7 = and ordered initial conditions with so = 10 and 
4> — 0.2. The time instances range from t = 2 10 to t — 2 17 . 
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FIG. 12. The size distributions of persistent empty regions 
(upper curves) and clusters (lower curves) for random initial 
conditions and 7 = 0.75. The measurement times are shown 
in the figure. 



E. Persistent Regions 

The mean-field arguments used for the empty site 
exponent predicted the region size distribution to de- 
pend on the initial cluster size distribution. Calcula- 
tions indicated a linear [see Eq. p|)] and an exponen- 
tial [Eq. (P9[)1 decay for the cases P(s ) — S(s — s ) and 
P(sq) = Sq exp(— s/so), respectively. It was further ar- 
gued that the dependence on the diffusion exponent en- 
ters only through a multiplicative factor of £~ 2 /( 2- t) m 
the distributions. This in turn means that the average 



region size would approach a constant at late times. Fig- 
ures [ll] and |l2] show the empty region distributions for 
these two cases and confirm the predictions. 

In Figure [ll] the smallest times show n are not large 
enough for the analysis of section [II C to be valid but 
the tendency of the distribution to approach a straight 
line is clearly visible. The reason why the asymptotic 
behavior is reached earlier in Figure [l^ is that the im- 
portant measure is not the time itself. Rather it is the 
ratio of the diffusion coefficients of clusters, which natu- 
rally grows faster the larger the diffusion exponent 7 is. 
Note, that in neither of the figures the horizontal axes 
have been scaled, which simply means that the average 
region size is constant, i.e., independent of time. This 
is consistent with the arguments presented in the sec- 
tion [V. The region size distribution between persistent 



filled sites behaves in the same way with a concentration 
dependent exponent 9f- 

As discussed at the end of section |^ the size distri- 
bution of persistent clusters remains unaltered both for 
monodisperse initial condition and for 7 = 0. Therefore 
we only present the result for random initial conditions 
and for 7 = 0.75 in figure [[2] The scaling is the same 
as it is for 7 = 0, i.e., the scaling function is a pure ex- 
ponential and the average cluster size is asymptotically 
independent of time. 



VI. CONCLUSIONS 

In this paper we have considered persistence in an ag- 
gregation process, in the simple case of one-dimensional 
DLCA. The emphasis is on local properties: cluster, 
empty site and filled site persistence together with the 
corresponding region and interval size distributions. We 
have shown that the three persistences are independent 
of each other with each behaving in its own, specific way. 

The perhaps most natural choice, the probability that 
a site has remained in the same state (either filled or 
empty), turns out to be nonuniversal. Nevertheless, only 
one of the two contributing probabilities, the filled site 
persistence, is responsible for this fact. The other one, 
the empty site persistence, is universal. In the DCLA 
the difference in the dynamics of empty and filled sites 
implies that there is no symmetry that would produce 
similarity in the persistence behaviors. The reason why 
the cluster persistence is a third independent quantity de- 
rives from the fact that it is a cluster property whereas 
the two other persistence definitions are considered at a 
fixed point in space. 

To summarize, the filled site persistence decays asymp- 
totically algebraically for 7 < 2 and for all concentra- 
tions: Pp(i) ~ t ■ The persistence exponent Of de- 
pends on concentration and is therefore nonuniversal. At 
low concentrations the filled site persistence decays as 
a stretched exponential upto a concentration dependent 
crossover time after which the collisions between clusters 



12 



become important and start to dominate the persistence 
behavior. 

The empty site persistence is universal. Although we 
have not been able to prove this, the universality is sup- 
ported both by mean-field type continuum arguments 
and by computer simulations. The former leads to a rel- 
atively simple relation between the persistence exponent 
and the dynamic exponent, 9e = 2z, verified by simu- 
lations. This is one of the few examples Q, where the 
inequality 9 > zd, where d is the spatial dimension, is 
fulfilled. The consequence of this is that the persistent 
empty regions do not have a fractal character. This is 
not true for example for the persistent regions in the 
Ising 0,0] or diffusion-annihilation jL7) models. The 
fact that 9e is notably larger than z, makes the separa- 
tion of the diffusive and persistence length scales clearly 
visible in the DLCA. 

The cluster persistence probability decays algebraically 
only for < 7 < 2 and as stretched exponentially for 
7 < 0. For < 7 < 2 the persistence exponent is 
given by 9c{l) — 2z = 2/(2 — 7) and it is discontinu- 
ous as 7 — > + since 9q (0) = |. All these results are in 
close connection with the scaling of the cluster size dis- 
tribution, especially with the small x decay of the scaling 
function f(x). In fact, there is a scaling relation between 
the exponents 9c = (2 — t)z. This independence of the 
cluster persistence probability on concentration and the 
strong support of the universal dynamic scaling behavior 
of the cluster size distribution J42[ speak for the univer- 
sality of the cluster persistence. The scaling relation to- 
gether with the result for 9c offers a new way to approach 
the determination of r for DLCA. 

It is worth emphasizing that the universal empty and 
cluster persistences decay with a same exponent 9e — 
9 C = 2z = 2/(2 - 7) for < 7 < 2 but have nothing 
to do with each other for 7 < 0. In fact, for 7 > we 
can write Pc(t) ~ [P£;(t)] r , where T takes the values 3/2 
and 1 for 7 = and < 7 < 2, respectively. A similar 
relation, P c (t) ~ [P E (t)] r , with a non-trivial T has been 
observed also for noninteracting random walkers [ fl9| and 
for the Potts model [H in one dimension. In DLCA the 
same persistence exponents, i.e., T = 1, for positive 7 
are due to the fact that the clusters giving the dominant 
contribution to the cluster persistence are those which re- 
main stationary. For other values of T the interpretation 
is more opaque. 

Although the persistences decay in a different man- 
ner (algebraic or faster) the intervals between persis- 
tent regions scale similarly for all cases: nx(k]t) = 
K x (ty 2 fx(t/K x (t)). This is due to the relation 9 X > 
z, which causes the persistent regions to be well sepa- 
rated at late times. Therefore the regions survive inde- 
pendently and the intervals size distribution is a simple 
exponential fx(y) — e~ y - 

The size distributions of persistent regions depend on 
the initial cluster size distribution. They do not depend 
on the diffusion exponent which only enters through a 
multiplicative time dependence on the distribution, for 



example, pE(l,t) ~ t~ 2 /( 2 ~^pE(l). This leads to a con- 
stant average region size at late times, which can be esti- 
mated from the knowledge of the initial cluster size dis- 
tribution. 

In conclusion, we have presented a rather comprehen- 
sive study of various local persistence probabilities in the 
one-dimensional DLCA. Our study is of interest also for 
the sake of practical realizations. Aggregation processes 
are plentiful, and all of the definitions - whether a point 
in space is occupied or not by a cluster, or whether clus- 
ters survive intact - might well be possible to measure in 
experimentally. It is interesting also to note that when 
the decay of a persistence probability turns out to be al- 
gebraic and universal, the exponent is always directly in 
some relation to the dynamical exponent z of the aggre- 
gation process. It is an obvious question to ask how the 
various quantities work out in higher dimensions. This 
generalization seems both challenging and interesting. 
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